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Abstract. An accurate description of the nuclear response functions for neutrino scattering 
in the Gev region is essential for the interpretation of present and future neutrino oscillation 
experiments. Due to the close similarity of electromagnetic and weak scattering processes, 
we will review the status of the scaling approach and of relativistic modeling for the inclusive 
electron scattering response functions in the quasielastic and A-resonance regions. In particular, 
recent studies have been focused on scaling violations and the degree to which these imply 
modifications of existing predictions for neutrino reactions. We will discuss sources and 
magnitude of such violations, emphasizing similarities and differences between electron and 
neutrino reactions. 



1. Introduction 

The motivation of the recent interest in neutrino-nucleus interactions is twofold. One one 
hand ongoing and future high precision neutrino oscillation experiments increasingly rely on 
our knowledge of neutrino-nucleus interactions in both Charged Current (CC) processes of the 
type 1^1 + A l + N+{A — l), where a virtual W boson is exchanged, and Neutral Current (NC) 
reactions I'l + A ^ i^'i + N + (A — 1), where a Z boson is exchanged. Typical neutrino energies 
involved in these experiments are of a few GeV. In this region a large contribution to the cross 
section comes from quasielastic and resonance processes, where the nuclear dynamics is known 
to play an important role. Since many electron scattering data exist at these kinematics (see 
Ref. m for a recent review), we should be able to use them to predict neutrino cross sections. 
As we will show, the superscaling properties of inclusive (e, e') data represent a bridge between 
electron and neutrino scattering reactions off nuclei. 

A second motivation consists in the fact that neutrinos can reveal informations on the nuclear 
structure complementary to electrons and they can be used to probe the internal structure of 
the nucleon (in particular the axial mass and the strange form factors), providing the nuclear 
content of the problem is under control. 

In this contribution we will first review the phenomenon of superscaling in inclusive electron- 
nucleus scattering O (Sj HJ [5] and then show some predictions for u-A cross sections based on 
the SuperScaling Approximation ("SuSA") [6l|7]. We will then discuss the microscopic origin 
of the superscaling function, with particular focus on two models, the Relativistic Mean Field 
(RMF) [8] and a BCS-like correlated Fermi gas [9], which seem to correctly reproduce the shape 
of the (e, e') experimental data, at variance with other relativistic models. Finally we will address 
the issue of scaling violations, in particular the ones connected to meson-exchange currents. 



2. Superscaling 

The experimental quasielastic (e, e') cross sections display the so-called "first kind scaling" [5] 
behaviour: at high momentum transfer {q larger than about 0.5 GeV/c) the reduced quasielastic 
cross section 

does not depend on two variables {q, co) but only on one combination of them, called the scaling 
variable and indicated by y{q, oj) in a non-relativistic scheme and ip{q, oj) in a relativistic context. 
These two scaling variables represent respectively the minimum momentum (— y) or kinetic 
energy (ttIj^V) that a nucleon inside the target nucleus must have in order to participate to the 
reaction. The conditions y = and 1^ = correspond, in the non-relativistic and relativistic case 
respectively, the quasielastic peak (QEP). In Ref. [4j it was shown that if the reduced function 
([1]) corresponding to different kinematical conditions is plotted as a function of y it fulfills with 
good accuracy first kind scaling in the so-called scaling region y < 0, whereas at y > scaling 
is broken. 

Second kind scaling consists instead in the observation that for different target nuclei the 
reduced cross section scales as the inverse of the Fermi momentum kp, so that the so-called 
superscaling function f = kp x F is independent of the specific nucleus. This has been illustrated 
in Ref. where it has been shown that scaling of second kind is excellent in the scaling region 
and slighthly violated at the right of the QEP. 

If first and second kind scaling are simultaneously respected the cross sections are said 
to "sup er scale" . A closer inspection of the Rosenbluth-separated longitudinal and transverse 
response functions shows that the scaling violations observed at ip > mainly reside in the 
transverse channel, essentially due to the A resonance contribution. 

Before proceeding, let us briefiy remind the basic formalism for inclusive electron scattering, 
focussing on the longitudinal channel where superscaling is seen to work better. In plane wave 
impulse approximation the longitudinal quasielastic response function can be expressed as 

RL{q,oj) = ^ / / dpd£S{p,£)nL{q,OJ-p,£), (2) 

q J Jt. 

where Ti{q,Lo) is the kinematically allowed region in the {£,p) plane, £ being the excitation 
energy of the residual nucleus and p its momentum, £) is the nuclear spectral function and 
TZl the single nucleon response, which depends on the electric and magnetic nucleon's form 
factors. 

Assuming that the single nucleon response TZl \s smoothly varying in the (f ,p)-plane and 
can be factored out of the integral one gets 

RL{q,uj) = TZLiq,uj;p = Pmm,£ = 0) X FLiq,uj) (3) 

where 

FLiq.Lo) = ^ / / dpd£S{p,£) (4) 

is the longitudinal scaling function. Note that the above factorization can be only approximately 
true in a relativistic context and it is violated to some degree due to off-shell effects and breaking 
of independent particle models. 

The simplest model where one can perform exact relativistic calculations is the Relativistic 
Fermi Gas (RFC), namely an independent-particle model in which on-shell nucleons are 
described by Dirac spinors. In spite of its simplicity the RFG model retains some fundamental 
properties: it is Lorentz covariant and gauge invariant and is therefore a good starting point 




Figure 1. The experimental superscaling function in the quasielastic (left panel) [TO] and A 
resonace (right panel) regions. 



for more sophisticated calculations. The RFG spectral function is S^^'^{p,£) = 9{kF —p)5{£ 
Tp + 7p) and the corresponding scaling function is 



q J Kkp J 4 

where the parabola f (1 — V'^) ^ (1 — V^^) = /(V') is the superscaling function and ^p = 'il:{q,uj) 
is the RFG scaling variable. The RFG model exactly superscales, as the data do in the so- 
called scaling region. However the comparison with the averaged data, illustrated in the left 
panel of Fig. [51 shows that the corresponding superscaling function is very different from the 
experimental one and in particular it fails the reproduce the long tail displayed by the data at 
large positive ip. However a simple phenomenological approach can be taken for describing the 
QE region: a fit of the experimental superscaling function can be extracted from the data and 
plugged into the RFG formulas in place of the parabola ([5]). Note that the fit involves only 
four parameters for describing all kinematics and all nuclei. Moreover, it represents a strong 
constraint on nuclear models, since any reasonable model used to describe the quasielastic region 
must reproduce this shape. 

The scaling analysis has been recently extended to the A resonance region [6]. In order to 
isolate the contributions in this region, the impulsive contributions arising from quasielastic eN 
scattering, calculated assuming superscaling is valid, has been removed from the experimental 
total cross section and the result has been divided by the elementary ^ A cross section, thus 
leading to a new superscaling function 

/a g,V'A) = kp r _j_ — T^aT • 

CTM [VLGf + vtG^I 

When the latter is plotted against the appropriate scaling variable V'A = ''PilP^'^p)-, which 
accounts for the different kinematics through the inelasticty paramater p = 1 + (m^ — 
m^)/(4|(5^|), the result shown in the right panel of Fig. [2]is obtained. Obviously this approach 
can work only at V'a < 0' since at > other resonances and the tail of DIS contribute. 

The two above superscaling functions can then be used to reconstruct the (e, e') cross sections 
in the QE and A regions, giving a satisfactory description of the data, as shown in Ref. [6]. 



3. Predictions for neutrino scattering: the SuSA approach 

We shall now illustrate how the superscaling properties of inclusive electron scattering can be 
exploited to predict neutrino cross sections off nuclei. 



Charged current neutrino-nucleus scattering can be easily linked to inclusive (e, e') reactions 
by replacing the virtual photon by a wirtual VF-boson. The main difference is that in the 
neutrino case, beyond the longitudinal and transverse response functions, a third response is 
involved, arising from the axial component of the weak current. In the RFG model the three 
response functions turn out to be _ _ 

Ri = G,x fiip) , (7) 

where the functions Gi depend on the weak nucleon form factors and the information about the 
nuclear dynamics is carried by the superscaling function /(V') • The basic assumption of the 
SuSA approximation is that the nuclear dynamics governing the neutrino-nucleus process is the 
same entering the electron scattering reactions and therefore the corresponding cross sections in 
the QE and A regions can be obtained by simply replacing in ^ the RFG superscaling function 
with the corresponding phenomenological fits: fRPG =^ fQE,A- 

In the left panel of Fig. [3] the double differential cross section in the QE and A regions is 
shown for (see Refs. [H] for more results): it clearly appears that in both regions the SuSA result 
lies significantly lower and extends over a wider range in k' than the RFG. In the right panel 
the SuSA fully integrated quasielastic cross section [11] on ^■^C is shown as a function of the 
neutrino energy and compared not only with the RFG but also with other relativistic models: 
the relativistic mean field (RMF) model, the relativistic plane wave impulse approximation 
(RPWIA) and a semi-relativistic shell model without (SRWS) or with (SRWS-tot) inclusion of 
discrete states of ^-^A^. Note that the RMF yields results which are very similar to the SuSA 
prediction, whereas the other models are closer to the RFG. The comparison of these results 
with the recent MiniBooNE data [12], which lie above the RFG and seem to point to a very 
high value of the axial mass {Mac^ = 1.35 GeV in place of the usual value 1.03, which we are 
using in our calculations), is somehow puzzling and still needs to be carefully explored. 




Figure 2. Left: CC neutrino reaction cross sections on C versus the outgoing muon 
momentum k', showing a comparison of the SuSA results (black online) with the RFG (red 
onleine); the separate QE and A contributions are shown with dashed lines. Right: integrated 
QE cross section versus the neutrino energy evaluated in different models. 

Concerning neutral current neutrino reactions, it has been proved in Ref. [7] that, although 
the kinematics is not the same of (e,e') reactions because the outgoing detected particle is the 
knocked nucleon instead of the lepton, the SuSA approach can still be used in the QE channel 
to predict cross sections. The result is qualitatively similar to the CC case, namely the SuSA 
predictions are lower the RFG and the response estends to a wider region. As anticipated an 
interesting feature of NC reactions is their sensitivity to the strangeness content of the nucleon, 
which allows in principle to extract informations on the strange hadronic form factors by different 
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Figure 3. NC neutrino (upper panels) and antineutrino (lower panels) cross sections for proton 
(left) and neutron (right) knockout. The effect of strangeness is shown in the various curves. 



combinations of neutrino and antineutrino scattering off protons and neutrons [13^ [TH [TS] . This 
is illustrated in the examples of Figs.[3l where it is seen that for proton knockout (left panels) the 
magnetic strangeness decreases (increases) the v (V) cross section, while the axial strangeness 
increases both v and 17 cross sections. For neutron knockout (right panels) and v scattering the 
situation is reversed with respect to p- knockout, while for V scattering the effect of magnetic 
strangeness is very small due to cancellation between the vector and axial-vector responses. 



4. Microscopic origin of the superscaling function 

In this section we will briefly illustrate the results of two relativistic models which, although 
sensibly different, seem to reproduce the basic properties of the experimental quasielastic 
superscaling function, namely a scaling behaviour of both first and second kind and a long 
high energy tail. 

The first model is the Relativistic Mean Field model, based on a Lagrangian containing strong 
scalar and vector potentials (corresponding to a, iv and p mesons) and describing the bound 
nucleon states as self-consistent Dirac-Hartree solutions. The final state interactions (FSI) are 
described in distorted wave impulse approximation using the same relativistic potentials [8]. As 
can be seen in Fig. H] (left panel) the RMF model provides the correct amount of asymmetry, 
unlike other relativistic models, also shown in the same plot. These correspond to the absence 
of FSI (RPWIA) and to a treatment of FSI through a real optical potential (rROP). It clearly 
emerges from the results that not only relativity, but also a consistent treatment of the initial 
and final states are required in order to get the observed asymmetry of the scaling function. 
Morevoer, as shown in Ref. [5], the model fulfills superscaling with good accuracy. 

A completely different model, which we will refer to as "BCS-like", has been shown to give 
very similar results [S] , as illustrated in Fig. [3] (right panel) . It is a simple covariant extension 
of RFG which accounts for NN correlations through the following wave functions 

1$ >= ]^(nfc + ?jfc4^aL^^)|0 > and [/^(p) >= ap| ]^[nfc + Ufc4^a^J [0 > (8) 
k l^pl k 

for the initial and final states, respectively. The Ufc-coefficients are taken phenomenologically 
as Fermi-Dirac distributions with one free parameter /3, controlling the tail of the momentum 




distribution v^. The associated superscaling function / turns out to be asymmetric in ip, the 
asymmetry increasing with the value of /?, and to reproduce the data for values of /3 ~ 0.02 
c/MeV. The superscaling properties of the BCS-like model have been explored in Ref. [9], where 
it was shown that both kinds of scaling are good at the QEP and broken to some degree at the 
right and left of it. 

5. Scaling Violations 

At low momentum and energy transfer (which can be defined approximately by g < 400 MeV/c 
and a; < 50 MeV) the scaling approach is bound to fail because collective eff^ects, like giant 
resonances, become important. However even at higher energies scaling violations are observed 
and it is very interesting to understand their origin. 

Among other possible contributions, like e.g. short range correlations, meson-exchange 
currents (MEC) can certainly play a significant role in breaking superscaling in the transverse 
channel. The MEC, being two-body currents, contribute to both the Ip-lh and the 2p-2h sectors 
of the response and their calculation is far from trivial in a relativistic context. Moreover in 
order to preserve gauge invariance, in parallel with the pure MEC diagrams, where the virtual 
photon attaches to the exchanged meson, it is necessary to take into account the corresponding 
correlation diagrams, where the photon attaches to the nucleon. Up to now an exact gauge 
invariant relativistic calculation has been performed only in the RFC system and in the Ip-lh 
channel |16 1I17[[T8] . Here the diagrams are the ones shown in Fig. O where thick lines in (d)-(g) 
represent either a nucleon or a A. It turns out that the diagrams involving the A give the 
dominat contribution and they break both kinds of scaling, as illustrated in the right panel of 
Fig. [5l Moreover a recent calculation [T9] performed in a semi-relativistic shell model essentially 
equivalent to the RMF with FSI has shown that the MEC contribution to the high energy tail 
of the response can become very large due to the dynamical nature of the pion propagator, an 
effect which is emphasized by the presence of FSI. 

In the 2p-2h sector the pure MEC have been calculated on the RFC basis [20j and shown 
to give a substantial contribution at high negative and positive values of ip, but a full gauge 
invariant calculation is still missing. 



6. Conclusions 

We have shown how the superscaling properties of (e, e') data can be exploited to predict charged 
and neutral current neutrino scattering cross sections and we have proposed two microscopic 
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Figure 5. Ip-lh MEC diagrams and results for / in the RFG (/o) and including the A-MEC 
diagram for different values of q (upper panel) and kp (lower panel). 

relativistic models, the RMF model and the BCS-like correlated Fermi gas, which are capable 
of reproducing the experimental superscaling function, in particular its long high energy tail. 

However, in order to assess the applicability of the SuSA approach and eventually improve the 
method we need not only to justify the microscopic origin of superscaling but also to understand 
its violations, in particular the ones associated with meson exchange currents. This has been 
partially accomplished in recent work [21j and more effort will be done in this direction in order 
to provide a consistent description of the nuclear dynamics in the quasielastic, resonance and 
deep inelastic regions. 
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